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ABSTRACT 


This  paper  derives  an  efficient  solution  procedure  for  solving  the 
continuous  version  of  the  Evans  modular  design  problem.  The  Kuhn-Tucker 
conditions  are  used  to  derive  a  dual  problem  which  can  be  solved  easily 
and  whose  dual  variables  indicate  which  equations  should  be  tight.  The 
technique  retains  a  tree-basic  solution  throughout  so  that  fast  solution 
routines  can  be  employed  which  are  quite  similar  to  tho:e  for  transporta¬ 
tion  problems.  Because  of  these  analogies,  the  solution  of  transportation 
size  problems  can  be  effected  with  only  moderately  increased  computer  time. 


1 .  INTRODUCTION 


The  modular  design  problem  was  first  presented  by  David  Evans  [4]. 
In  this  problem,  parts  are  to  be  grouped  into  a  single  module,  several  of 
which  are  then  used  in  each  application.  The  objective  is  to  minimize  the 
total  cost  of  parts  used: 


where 


Min  E  c,  x'  £  d,  y' 
*,y  lei  1  1  jcj  -1  J 


S  .  t. 


/  /  ^  / 
x.  y .  >  r .  . 
i  J  ~  ij 


integers 


|  For  all  i  and  j 

r 


i 

j 

c . 

1. 

d 

j 


t 

X. 

1 


/ 


[l >  2 , . . . ,m] 

[l , 2 , . . . , n  ; 

cost  of  part  i 

demand  for  application  j 

number  of  part  i  units  required  in  application  j. 

the  number  of  part  i  on  the  module 
(decision  variable) 

the  number  of  modules  needed  in  application  j 
(decision  variable) 


In  the  previous  approaches  to  this  problem  listed  in  the  bibliography 
and  also  in  the  present  paper,  the  problem  is  modified  to  the  continuous 
modular  design  problem  by  dropping  the  integer  requirements  on  and  y'  for  all 
i  and  j;  it  is  believed  that  a  solution  to  the  continuous  problem  is  a  necessary 
prelude  to  any  solution  of  the  integer  version. 
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After  the  Evans  paper,  a  second  paper  on  modular  design  was  written 
by  A.  Charncs  and  M.  Kirby  [2].  Both  of  these  papers,  proposed  solution 
procedures  based  on  searching  the  x-y  space  via  specialized  search  routines. 

A  third  paper  l 6],  by  A.  Passy,  modified  Charnes  and  Kirby's  procedure  by 
formulating  the  model  as  a  geometric  programming  problem.  The  approach 
p’-o^^nted  in  Persy's  paper  is  tr  ~cvc  from  one  system  of  tight  constraints 
to  another  until  the  optimum  is  found.  Although  Passy's  procedure  avoids 
the  relatively  slow  search  procedures  of  the  first  two  papers,  it  has  three 
major  drawbacks.  First,  convergence  of  the  procedure  was  not  proved.  (The 
results  of  the  present  paper  might  aid  in  proving  convergence  for  Passy's 
procedure.)  Second,  procedures  outlined  by  Passy  to  solve  the  problem  of 
cycling  could  (and  most  likely,  would)  lead  to  an  inordinate  number  oi  pivots. 
And,  finally,  the  complexity  of  '"he  calculations  for  each  pivot  step  could 
require  excessive  amounts  of  computer  time  even  if  convergence  were  proved. 

It  is  the  intent  of  the  present  paper  to: 

1. )  Develop  a  dual  problem  (with  properties  similar  to  the  linear 

programming  dual)  from  the  Kuhn-Tucker  optimality  conditions. 

2. )  Use  the  dual  solutions  to  develop  a  simplex-like  solution  algorithm 

and  prove  its  convergence. 

Preliminary  tests  indicate  that  the  algorithm  presented  in  this  paper  is 
very  efficient.  The  authors  have  found  that  a  continuous  modular  design  problem 
can  be  solved  with  only  slightly  more  computation  time  than  a  transportation 
problem  of  the  same  size.  They  are  presently  working  on  further  computational 
tests  of  the  procedure  which  will  be  reported  on  elsewhere. 

2.  DUALITY  IN  THE  MODULAR  DESIGN  PROBLEM 

Kvans  modified  the  '’riginal  problem  by  making  the  following  substitutions: 
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y  •  =  d  4  y  ' 
j  J  j 


r.  .  =  r. .  c.  d. 
iJ  ij  i  j 


He  also  noted  that  there  exists  an  infinite  number  of  solutions  to  the  problem 

since  if  x  and  y  are  solution  vectors,  then  so  are  x/6  and  y»0  for  any 

3  >  0.  We  may,  therefore,  add  the  restriction  >’y.  *1  which  singles  out 

j  2 

a  unique  member  of  the  class  without  loss  of  generality.  These  transformations 
lead  to  the  primal  problem  (P): 

(IP) 


min  I  x.  =■  g 
x  i 


s.t.  x.y .  -  z.  .  «  r. .  for  all  i  and  1 
i  J  ij  ij  J 


sy.-i 
j  J 


z.j.  ,  y j ,  r_  >  0  For  all  i  and  j 
where  z.  .  is  a  surplus  variable. 


(2P) 

(3P) 

(4P) 


Theorem  1  (Evans)  The  solution  to  problem  (P)  exists  and  is  unique. 
This  theorem  was  first  proved  by  Evans  {4] .  The  algorithm  established 
in  the  present  paper  gives  an  alternate,  constructive  proof  of  the  existence  of 
a  solution  to  problem  (P). 

LEMMA  1.  If  the  requirements  matrix,  R,  has  no  zero  revs  or  columns, 

then  in  any  solution,  x^,  y^  >  0  for  all  i  and  j. 

Proof.  The  assertion  is  obvious  since  x.y,  >  0  at  least  once  for 

i  j 

each  i  and  j. 

Let  X  be  associated  with  constraints  (2P)  and  p,  bo  associated  with 
constraint  (3P).  Then  it  is  easy  to  show  that  the  Kuhn-Tucker  conditions 
associated  with  the  primal  problem  are: 


-4- 


f  Vi 

*  1 

for  all  iel 

(1C) 

fVu 

*  P- 

for  all  jej 

(2C) 

z .  .  •  X ,  . 
ij  iJ 

=  0 

For  all  i  and  j 

(3C) 

X.  .  >  0 

ij  - 

For  all  i  and  j 

(4C) 

LEMMA  2.  For  any  pair  of  feasible  solutions  to  the  primal  and  dual 
problems,  p,  *  g. 

PROOF.  Multiply  (1C)  by  x^ ,  sum  over  all  i,  and  use  (IP)  to  show: 

£  E  xi  X  y  -  £  xt  -  g 

i  i  J  i 

Now  multiply  (2C)  by  y^,  sum  over  all  j  and  use  (3P)  to  show: 

E  E  Xf  X  y  *  U  ly,  =  u 

i  j  J  1 


Hence,  p  *  g. 

By  analogy  with  classical  linear  programming  we  shall  interpret  con- 
striints  (3C)  as  "complementary  slackness"  conditions  and  insure  that  they 
hold  by  the  algorithmic  solution  techniques  we  develop.  Again  by  analogy  to 
linear  programming  we  add  the  objective  function  E  E  X^ ,  r^ .  to  constraints 


(1C),  (2C) ,  and  (4C)  to  create  the  dual  problem.  The  remainder  of  this  section 
is  devoted  to  showing  that  most  of  the  simplex  method  solution  techniques  involv¬ 
ing  interplay  between  the  primal  and  dual  problems  can  be  carried  over  to  the 
modular  design  problem  considered  here.  In  later  sections  we  show  that  they 
are  powerful  enough  to  make  possible  an  efficient  simplex-like  algorithm  for 
solving  the  continuous  modular  design  problem. 


The  dual  problem  is  defined  by: 


Maximize 

E 

V 

X.  ,  r 

..  -  f 

(ID) 

X.  . 

1 J 

i 

J 

ij 

ij 

t:'U  'J 

- 

L 

tui  all  i 

(2D) 

s.  t. 
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I  x.  X^ .  ■  g  for  all  j  (3D) 

i  3 

X^  >  0  for  all  i  and  j  (4D) 

LEMMA  3.  For  any  pair  of  primal  dual  solutions  (whether  nonnegative 
or  not) :  g-f=‘£Z\..*z,1. 

,  .  Lj  lj 

1  j  J 

PROOF.  Multiply  (2P)  by  sum  over  all  i,j  and  use  (ID)  to  show: 

1 1  *i '  f  Vli  r‘J‘£ 

but  S  E  L .  x  y.  *  g  (see  proof  of  Lemma  2) 

i  j  1J  3 

hence,  g  -  f  3  EEL.  z  .  . 

i  j  3  3 

LEMMA  4.  (Complementary  Slackness.)  For  any  pair  cf  primal  dual  feasible 
solutions  g  *  f  iff  X^-z^  =  ®  *or  flXl  i,j. 

PROOF.  The  proof  follows  directly  from  Lemma  3  and  the  fact  that  X^  and 
are  non-negative  for  all  i  and  j. 

THEOREM  2  (Duality  Theorem).  The  quantities  ,  y^  for  all  i  and  j 
are  a  solution  to  problem  (P)  iff  X^  ^  for  ail  i  and  j  are  a  solution  to 
problem  (D)  and  g  =  f. 

PROOF.  (The  proof  of  this  theorem  will  be  only  briefly  sketched.] 

The  Arrow-Hurwicz-Uzawa  Constraint  Qualification  (5.  p,  102]  can  be  shown  to 
hold  for  the  constraint  set  of  the  primal  problem.  This  implies  that,  at  the 
optimum  solution  to  the  primal  problem  there  exists  a  solution  to  the  Kuhn-Tucker 
problem  (5,  pp.  105-106] .  The  functions  defining  the  constraint  set  can  be 
shown  to  be  quasi -concave  in  the  non-negative  orthant  while  the 
objective  function  is  linear.  Zangwill  [12 ,  p.  43]  has  shown  that  in 


such  a  case  a  solution  to  the  Kuhn-Tucker  problem  can  occur  only  at 
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the  optimum  of  the  primal  problem.  Finally  the  K-T  problem  is  made  up 

of  the  constraint  set  of  the  dual  problem  together  with  the  complementary 

slackness  conditions  *z.,  *>  0  which  are  true  if  and  only  if  g  =  f 

ij  ij 

(Lemma  4),  This  completes  an  outline  for  the  proof  of  Theorem  2. 

A  property  of  interest,  although  not  used  here  is  that  problems  (P) 
and  (D)  are  not  mutually -dual.  If  we  form  the  dual  oL  problem  (D)  we  get 
the  following  problem: 

g'  (IE) 


Min 

^>3 


iqnn. 
i  j 


s.t.  x.  T(.  -e  C.y.  >  r..  for  all  i  and  \ 
l  j  *Vj  -  ij 


(2E) 


where,  C.  is  the  dual  variable  associated  with  constraint  (2D)  and  Tt  is 

’  j 

the  dual  variable  associated  with  constraint  (3D). 

For  fixed  x.  and  y.  problems  (D)  and  (E)  are  mutually  dual 
-  J 

generalized  transportation  problems.  If  and  y  are  solutions  to  (P) 

then  s  xi  ’  =  ®  are  a  feasible  solution  to  (F.)  from  which  it  easily 

follows  that  g  >  g 

3.  THE  RESTRICTED  PROBLEM  AND  ITS  SOLUTION 

In  this  section  it  is  first  necessary  to  define  some  properties  of 
graphs.  For  a  general  discussion  of  graph  theory  see  Berge  [lj.  The  follow¬ 
ing  description  was  taken  from  [10,  p.-2], 

"Let  V  be  a  set  of  n  elements  called  vertices  or  nodes  and  let  E  be 
a  set  of  (some  of  the)  pairs  (u,v)  with  u,v  eV.  A  pair  (u,v)  is  called  an  edge 
between  u  and  v,  or  also  between  v  and  j  (no  direction  is  Implied). 

Then  G  =  (V,E)  is  called  a  graph .  A  path  between  u  and  v  in  G  is  a  list 
U  *  W  w  . w  =  v 

»'  t  I 


where,  (w^ wj)  e  E  for  J  =  1.  •••»*:.  A  path  is  a  cycle  if  u  =  v  in  the 
above  list.  A  graph  is  acyclic  if  it  has  no  cycles.  A  graph  is  connected  if 
there  is  at  least  one  path  connecting  each  pair  of  distinct  nodes.  A  tree 
is  a  connected  acyclic  graph.  Equivalently,  a  graph  is  a  tree  if  and  only  if 
there  is  a  unique  path  between  each  pair  of  distinct  nodes." 

In  addition  to  the  definitions  quoted  above  we  shall  need  the  following. 
A  forest  is  an  acyclic  graph.  It  is  easy  to  show  that  a  forest  is  the  union  of 
trees,  that  is,  a  union  of  connected  acyclic  graphs. 

In  the  modular  design  problem  we  shall  consider  the  graph  G  *  (V,E) 
defined  as  follows:  The  set  V  of  nodes  consists  of  the  rows  and  columns 
of  the  requirements  matrix  R;  the  set  E  of  edges  consists  of  some  of  the 
cells  (i,j)  of  the  R  matrix. 

Suppose  C  has  a  cycle 

^  —  » ^2  ^  ’  *  *  *  *  ^ 

where,  sp  =  sp+1  or  t?  =  t  ^  for  p  = 


(Note  ■  Sj  and  t^+^  *  t^)  and  L  is  an  even  number.  In  each  row  or  column 

of  the  R  matrix  there  are  either  zero  or  two  cells  of  the  cycle.  Then  I'  can 
be  written.  I'  =  U  where: 


ri  *  ((s^tj),  (s.^.tj), - (s/v.x> 

r2  =  {(s2,t2>,  (s4,t4),...,(s^,t^)} 


DEFINITION.  The  value  w, 

-  (s,t) 


Is  defined  to  be  the  ratio 


of  a  cycle  relative  to  any  element  (s,t)  c 


DEFINITION.  A  nondegenerate  problem  is  one  having  w 


*  1  for 


(s,t) 


all  cycles  f  and  (s,t)  e  Tj. 


LEMMA  5.  A  problem  with  data  r^ .  for  all  i  and  j,  may  be  replaced 

by  a  nondegenerate  problem  with  perturbed  data  +  6*+m^  if 

r..  0  or  r.  =  0  if  r.  .  =  0  and  where  5  can  be  chosen  arbitrarily 

iJ  ij 

sm  ill. 

PROOF.  For  given  m  and  n  there  are  only  a  finite  number  of  possible 
cycles  T.  For  such  a  cycle  to  have  value  1  with  the  perturbed  data  we  must 
have 

i  1  in  4 

r  (r 


(r . .  +  6i4TCj) 


<  i ,  j  )  e  i" 


1 


(s,t)«r. 


'  st 


ss4tnt)  =  0. 


This  expression  is  a  polynomial  in  b.  Moveover  there  is  at  least  one  power 
of  C  that  has  a  non  zero  coefficient.  To  show  this,  let  (h,k)  he  the  cell  in 
I  with  smallest  k,  and  given  this  k  the  smallest  h;  suppose  ( 1 1 , k ) s (a 
similar  proof  holds  for  (h.k)?!'^).  Then  there  is  a  term  c  t  where 
c  =  i .  #  0  since  the  value  of  the  cycle  is  one  and  all  other 

fi,j>er,-{(h,k)}  11 

tf  rms  have  higher  powers  of  b.  Hence  we  need  only  choose.  5  small  and  not  equal 
t.o  any  of  a  finite  number  of  zeros  of  a  finite  number  of  polynomials  to  obtain  a 
nondegeaerate  perturbed  problem  close  to  the  original  one. 

In  the  rest  of  this  paper,  we  shall  assume  that  we  are  dealing  with  a 
nondegenerate  problem.  The  techniques  for  extending  the  algorithm  we  shall 
present  to  degenerate  problems  are  similar  to  those  for  linear  programming 
and  will  not  be  discussed. 

DEFINITION.  Given  a  feasible  solution  to  problem  (P)  by  a  tight  constraint 
in  row  i  we  shall  mean  a  cell  f i  .  j )  such  that  "  rjj  “  Zj  -j 
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LEMMA  6.  At  an  optimum  solution  of  problem  (P)  there  is  at  least  one 

tight  constraint  in  each  row  and  column  of  R. 

PROOF.  Assume  the  contrary,  that  we  have  an  optimal  solution  and  for 

some  row  u  no  cell  is  tight,  i.e.,  x  y.  >  r  .  for  all  j.  But  then  we 

u  j  u  j 

can  decrease  x^  while  keeping  the  solution  feasible  and,  therefore,  decrease 
the  objective  function,  which  is  a  contradiction.  By  the  problem  symmetry, 
the  same  kind  of  proof  is  valid  for  columns  of  R. 

LEMMA  7.  If  x. ,y j ,  and  ^  are  solutions  to  problem  (P)  then  the 
graph  with  nodes  being  rows  and  columns  of  R  and  edges  being  the  tight  cells 
(i.e.,  z  =  0)  is  a  forest. 

PROOF.  The  proof  follows  directly  from  Lemma  6  and  the  fact  that  we 
are  dealing  with  only  non-degenerate  problems,  therefore  eliminating  the 
possibility  of  cycles. 

We  shall  now  define  a  forest- res  trie ted  problem  associated  with  a  given 


S  .  t. 


to  be : 

rain  Z  x, 

Vyj  i€l  1 

=  % 

(IH) 

i  yj  -•  1 

jeJ  J 

(2H) 

Vj  -  rty 

(i,J>  s  F 

(3H) 

Vy j  *  0 

for  all  i  and  j 

(*H) 

A  tree  restricted  problem  is  a  forest-restricted  problem  where  the  forest  is 
made  up  of  a  single  tree.  In  the  algorithm  to  be  presented  in  this  paper 
the  authors  retain  a  tree  basis  throughout  the  procedure.  We  are.  therefore, 
primarily  interested  in  solutions  to  tree-restricted  problems.  The  solution 
to  a  tree-restricted  problem  will  now  be  characterized. 

For  any  tree-restricted  solution  there  exists  a  unique  path  connecting 
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any  two  columns  of  R  f l] ,  f 10] .  This  unique  path  may  be  presented  as  in 
Figure  1,  for  the  case  that  column  1  is  connected  to  column  q.  We  have 
displayed  only  those  columns  aid  rows  from  R  that  correspond  to  the  path 
between  column  q  and  column  1.  In  Figure  1  the  cells  where  r's  appear 
are  all  tight  and  rows  and  columns  have  been  permuted  and  relabeled  to  be  in 
the  staircase  form  showed.  Certain  other  cells  may  be  tight  in  these  rows  and 
columns  but  are  not  of  interest  at  this  time,  hence  they  are  not  indicated  in 
Figure  1. 


rll 

r12 

r22 

r23 

• 

•  r  i  i 

p-l,q-l 

1 

H 

r 

pq 

i 

! 

. 

Figure  1 

It  is  easy  to  see  that 


y 


q 


r 

_ 2-0 

r  . 

p»q-i 


q-i 


since  x  y 

p  q 


r 

pq 


and  x  y  ,  *  r  The  procedure  can  be  continued  in  a  similar  fashion  until 

P  q~ 1  P,q-l 
finally : 

y  =  _l£s _  •  ...  In.  ■  In 

q  rp.q-l  rp-l,q-2  r22  rll 


yi "  Vi  (1K) 
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In  general  denote  by  d  the  ratio  ;_u  •  It  WH1  always  be  the 

yt 

case,  as  above,  that  d^  is  the  quotient  of  products  of  r^'a.  We  nan 

choose  an  arbitrary  column,  say  column  k,  in  the  matrix  R  and  represent 

all  values  of  y  for  jeJ,  in  terms  of  by  means  of  the  equation 

y.  =  y^.  Note  th  it  d^j.  *  1.  Using  constraint  (2H)  we  have  the  solution  to 

the  tree  restricted  problem  as  I  y.  =  Z  d ..  *y.  ■  1  which  yields 

jeJ  J  jsJ  J 


I  d.. 

Jk 


and  y  =  d  ,y,  for  all  veJ.  The  values  of  x,  are  easily  determined  from 
Jv  vk  k  i 

constraint  (3H). 

Given  a  tree  basis  T,  the  associated  restricted  tree  dual  solution 
can  be  derived  by  using  the  fact  that  X . .  *  0  for  (i,j)  not  contained  in 
the  tree  basis.  This  follows  since  these  cells  are  not  forced  to  be  tight. 
Also,  .  =  0  for  (i,j)  contained  in  the  tree  basis.  Taking  constraint 
(2P),  multiplying  by  X  /g  and  summing  over  i  and  j  separately  yields: 


£  X.  .  x. 
i  1 

y./E  -  tXij-ij/s- 

7  X  .  ,  *r .  ,/g 

L  i.)  ij  & 

and 

Z  X.  .  x. 
i  iJ  1 

yi/s  '  ‘hj'V8  ’ 

Z  X . . *r  /g 
j  i]  i.l 

using  constraints  (2D) 

and  (3D)  and  complementary  slackness 

we  may 

equations  to 

he : 

yj = 

<i,j)«T  ij  11  8 

for  j  e J 

(11.) 

and 

xi  * 

(i,j)eT 

for  iel 

(2L> 

It  is  interesting  to  note  that  the  form  of  the  equations  for  the  variables  X. 


ij 


are  the  same  as  the  form  derived  by  Pussy  [6,  p.  450]  for  the  geometric  dual 


* 
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variables  if  we  let  =  \^*r„/g.  The  dual  objective  function  developed 

by  Passy,  however,  is  very  different  from  the  one  used  in  this  paper. 

The  tree  structure  of  the  nonzero  variables  in  equations  (1L)  and  (2L)  means 
chat  the  equations  can  be  solved  by  a  simple  solution  procedure.  Kote  that  they 
can  be  rewritten  as 


yj = 

r.  Pij 

(iJ)eT  1J 

for  jeJ 

(ll') 

/ 

xi  - 

S  P^ » 

<i,j)eT 

for  iel 

(2L') 

where  we  have  made  the  substitutions  x.'  *  x./g  and  p.  .  “  X,  .  r.  ,/g.  The 

r  r  ^i3  ij  ij 

following  algorithm  finds  the  p_'s  given  the  primal  solution  x^.y.: 

(1)  Let  TR(TC)  be  the  3et  of  rows  (columns)  containing  a  unique 
tight  cell.  Because  we  have  a  tree-restricted  solution  TRUTC  is  not  empty. 

(2)  For  all  tight  cells  (i,j)  with  ieTR  let  p.^  *  x^.  For  all 

tight  cells  (i,j)  with  jeTC  let  *  y^ .  That  this  is  correct  follows 

from  the  fact  that  these  tight  cells  are  unique  in  their  rows  or  columns. 

(3)  "Cross  out"  the  rows  i  eTR  and  j  sTC.  For  the  remaining 
matrix  define  a  new  primal  solution  x^  and  y^  as  follows 

x  '  «=  x /  -  I  y.  for  all  i  itR 
jsTC  J 
U,j)€T 

y.  *  y.  -  Z  x'  for  all  j  ^TC 
J  J  ieTR 
( i  , j ) eT 

(4)  if  there  are  no  uncrossed  out  rows  stop;  otherwise  go  back  to 
step  (1)  and  repeat. 

A  few  additional  properties  of  tree-restricted  solution  will  now  be 
presented . 
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LEMMA  8.  For  any  tree  restricted  solution  let  X  be  the  tree  basis 
and  let  (p,h)  be  a  ceil  not  in  T.  Then  Tll{(p,h)}  has  a  unique  cycle 
whose  value  is 


w. 


x  yu 
-P-Jl 


(P,h>  r?h 

PROOF.  Let  us  assume  (without  loss  of  generality)  that  h  *  1  and 
column  1  is  connected  to  row  p  as  in  Figure  1,  Adding  the  cell  (p,h)  *  (p,l) 
we  obtain  the  cycle  shown  in  Figure  2  .  We  know  that  y^  *  d^y^  and 


1 

2 

•  « 

•  % 

0 

1 

rll 

r12 

•  •  *  • 

»  »  •  • 

2 

• 

r 

r22 

•  «  *  • 

« 

rp-l,q-l 

P 

r  , 

p,q-i 

r 

p.q 

« 

_ 

* 

* 

. 

• 

_ 

[ _ ! _ 1 

Figure  2 


also  that  x  =  r  /y  .  Hence 

p  pq  q 


Vi 


r 

_E2 

V 


where,  as  in  (IK), 


Jjl zhszi 

rp,q-l  rp-l,q-2 


xa,. 


r23 

C 12 

r22 

ril 

i' 


It  follows  that 
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x  y,  r  r  .  r 

p  i  _  _&a  jLxSzl  _r 


r22  rll 


r  i  r 

pi  pq 


p-l,q-l 


(p,l)* 


as  wa«  to  be  shown. 

COROLLARY.  A  tree-restricted  solution  x^>Yj  primal  feasible 
if  and  only  if  every  non  basic  cell  determines  a  cycle  whose  value  is  >  1. 

This  corollary  is  used  to  check  for  primal  feasibility  in  the  algorithm 
to  be  presented. 

DEFINITION.  Given  a  tree  basis  T  and  any  cell  (p,q)eT  we  define 
the  following  four  sets: 

I  =  [set  of  all  rows  that  can  be  reached  in  T  using  cell  (p,q), 
except  for  row  p{ 

I  =1-1 
P  q 

Jp  =  {set  of  all  columns  that  can  be  reached  using  cell  (p,q), 
except  for  column  q} 

J  =  J  -  J 

q  P 

Clearly  pel  and  q  e  J  and  those  sets  are  never  empty.  Also  at  most 

P  q 

one  of  the  sets  I  and  J  is  empty, 

q  P 


J  J 

p  q 


•<p,q)  : 

Figure  3 


Figure  3  shows  the  matrix  R  divided  into  the  four  subsets  1  x  J  , 

P  P 

I  x  J  ,  I  x  J  ,  and  I  x  J  .  Of  these  four  sets;  I  x  J  contains  no 
pqqp’qq  qp 

cells  of  T;  1  x  .J^  contains  only  the  cell  (p,q)  €  T;  all  the  rest  of 

the  cells  of  T  are  in  the  other  two  areas  (I  x  J  )  U (I  x  J  ). 

P  P  q  q 
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In  the  algorithm  to  be  presented  later  we  are  going  to  change  data 
elements  r  in  a  parametric  fashion.  The  next  two  lemmas  characterize 

pq 

what  happens. 

LEMMA  9.  Let  x. ,  y ^  be  a  primal  feasible  tree-restricted  solution 

with  tree  basis  T  and  let  (p,q)eT.  If  we  replace  r  by  a  larger  value 

r  >  r  ,  then,  provided 
pq  -  pq 

ic 

r  <  Min  r  w 

pq  -(=,0SIqxJp  ”  (s'c) 

the  tree-restricted  solution  x^ ,  with  the  same  basis,  T,  is  primal  feasible. 

PROOF.  By  the  corollary  to  Lemma  8  we  need  only  show  that  every  non 

basic  cell  (u,v)  determines  a  cycle  with  value  >  1.  Referring  to 

Figure  3  it  is  obvious  that  non  basic  cells  in  the  (I  x  Jp)  or  (1^  x  J^) 

areas  have  cycles  entirely  contained  in  these  areas.  Hence  changing  r 

does  not  affect  the  values  of  their  cycles. 

Because  every  cycle  goes  alternately  from  row  nodes  to  column  nodes,  every 

non-basic  cell  (u,v)  in  the  Ip  x  area  determines  a  cycle  F  which  includes  (p,q) 

in  the  Tj  and  (u,v)  in  the  I'2  part.  Hence  w(y  v)  increases  if  rp<j  increases  and 

* 

primal  feasibility  continues  to  hold  for  these  cells  however  large  we  make  rp^. 

Finally  consider  cells  (s,t)  in  the  I  x  Jp  area.  Such  cells  determine  a 

cycle  F  with  (s,t)  and  (p,q)  in  the  7^  part.  Let  and  w^g  ^  ke  value 

of  the  cycle  determined  by  ts,t)  with  r  and  rp^  >  rp^  respectively. 

Then  w*  ,  r  "  =  w,  .  r  by  the  definition  oi  the  value  of  a  cycle. 

(s,t)  pq  (s,t)  pq 

Since  wc  want  w.  ,  >  1  this  means  that  we  must  have  r  <  w  ,  r  tor 
(s  ,  t )  -  pq  Vs »  l  ;  pq 

every  (s.t)  e  l  x  J  and  therefore  the  statement  of  the  Lemma  is  true. 

J  q  p 

LEMMA  10.  Let  x^,  be  a  primal  feasible  tree-restricted  solution 

* 

with  the  tree  basis  T  and  let  (p,q)eT,  If  we  replace  rp^  by  n  smaller 
**  * 

value  r  r  , 
pq  -  pq’ 


then  provided 
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r  * 

**  .  Ipq 

r  >  Max 

pq  (u.vJelpXJ^-CCp.q)}  W(u,v) 

the  tree-restricted  solution  ,  y^  with  the  same  basis,  T,  is  primal  feasible. 

PROOF.  The  proof  here  is  analogous  to  that  for  Lemma  9  and  will  not 
be  given. 


4.  THE  MODIFIED  PROBLEM 


Given  problem  R  with  data  r„  we  define  the  modified  problem  R 
★  * 

with  data  r.  .  where  r  .  <  r  <  ®  for  all  i  and  j.  We  shall  give  an 
*■  J  ^  3  *•  J 

* 

algorithm  for  finding  an  optimal  tree-restricted  solution  to  R  and  show 
how  this  can  be  used  to  find  the  optimal  forest  basis  solution  to  Problem  R. 

it 

THEOREM  3.  Given  an  optimal  tree-restricted  solution  to  R  there 

corresponds  a  unique  forest-restricted  optimal  solution  to  R.  Conversely, 

to  the  optimal  forest-restricted  solution  to  problem  R  there  correspond  at 

* 

least  one  optimal  tree-restricted  solution  to  R  . 

*  *  * 

PROOF.  Given  an  optimal  tree-restricted  solution  x.  ,  y.,  X  .  ,  with 

i  J  ij 

tree  basis  T  to  problem  R  drop  from  T  all  tight  cells  such  that  X^  =  0. 

ic  Vc 

The  result  is  a  forest  F  with  X.  .  >  0  for  (i,j)  e  F  and  X.  .  11  0  for  (i,i)  t 

ij  ij 

*  'k  ★ 

Since  >  r^  >  r^  for  all  i  and  j  these  solutions  are  primal  and 

dual  feasible,  and  hence  by  the  duality  theorem  are  optimal  for  problem  R. 

Given  an  optimal  forest-restricted  solution  x. ,  y.  and  X..  with  forest 

1  J  ij 

basis  F  we  shall  give  a  constructive  procedure  for  deriving  an  optimal  forest 

* 


restricted  solution  to  a  problem  R  .  Suppose  F  =  U  U  ...  li  T^  where 
each  T.  is  a  tree.  If  row  u  contains  a  (tight)  cell  of  T then  it  will 
not  contain  a  cell  from  any  other  tree.  Similarly,  if  column  v  contains  a 
cell  of  T  then  it  will  not  contain  a  cell  of  any  other  tree.  We  now  show 
one  way  to  "hook  together"  the  trees  in  F  and  make  them  into  a  single  tree. 


F. 
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Reproduced  from 
best  available  copy. 


Let  be  the  index  of  any  row  containing  a  cell  of  and  let  »Jr 

be  indices  of  columns  containing  cells  of  Add  the  cells 

^1*^3)  ’  ‘  *  *  *  ^l’^fc)  t0  F  whicfl  ^11  make  it  into  a  connected  tree  aod  also 
*■ 

define  problem  R  by 


W 

y  4  -  ri  1  *  *  *  * ,ri  i 

L1  J2  ,J2 


•ff 

and  all  other  r  ■  rij*  ^°^ows  t^iat  x£>yj 

★ 

dual  feasible  and  hence  optimal  for  R  . 


and  X 


ij 


ate  still  primal  and 


Obviously,  there  are  many  other  ways  the  trees  . . . tT^  may  be 

connected  to  male  a  single  tree  so  the  above  process  is  not  unique. 

* 

THEOREM  4.  Let  x^,  y^  be  a  feasible  solution  to  problem  R  with 
tree  basis  T  and  let  be  the  restricted  dual  solution;  and  consider 


a  cell  (p,q) 

e 

T ;  then 

(a) 

if 

X  <0 

pq 

we  can  decrease 

* 

g  by  increasing 

* 

r 

pq 

(b) 

if 

X  >0 

pq 

we  can  decrease 

■k 

g  by  decreasing 

* 

r  . 

pq 

These  results  hold  only  over  a  sufficiently  small  range. 

PROOF.  Suppose  we  set  r  *  r^  +  6^  and  write  fhe  Lagrangian 
function  of  the  primal  problem.  It  is 


L(*,y,X,n,z, 6) 


I  E  y.  +  -  (rtJ  +  61j)I 


+  n(  E  y,-l) 

J  J 

Now  holding  all  6^  v  fixed  at  zero  except  for  6  and  letting  g(5jj)  be  the 
corresponding  value  of  the  primal  problem  we  can  rewrite  the  Lagrangiajl(for 
small  changes  in  6^)  ae 

8(&ij)  *  8(0)  +  Xjj  +  0(6^) 


for  which  the  two  assertions  are  obvious. 


1 
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ic  k  it 

THEOREM  5.  Giveu  a  feasible  tree-restricted  solution  x^, 
to  problem  R  with  basis  T,  let  (p,q)  s  T.  Define  r°  to  be  the  value  of 

r  t 

k  k 

r  for  which  >  =  0,  then 

pq  pq 


where  the  quantities  e^p,  d^,  Rp»  and  will  be  explained  in  the  proof 
below.  Then 


(a) 

X 

* 

o 

X  = 

G 

«  r  = 

r 

pq 

pq 

pq 

(b) 

★ 

0 

* 

o 

x  < 

«  r  < 

r 

pq 

pq 

pq 

<«> 

* 

* 

o 

X  -' 

0 

«=>  r  > 

r 

pq 

pq 

pq 

PROOF 

Let 

the  sets 

i 

‘"'jaihu iZ,,. 


p'  q  p  q 

we  remove  (p,q)  from  T  then  as  noted  before  all  the  remaining  cells  of  T 


are  in  the  areas  I  x  .1  and  I  x  J  . 

P  P  q  q 

For  each  j  e  Jp  let  i  be  the  smallest  row  index  such  that  (i,j)  e  T; 

* 

and  let  Rp  be  the  set  of  such  cells  (i,j).  Then  ■  riyxi  for  (i,J)  e  Pp 


Also  let  e,  be  the  ratio  e, 
iP  ip 


x.  /x  . 
i  P 


Similarly  for  aach  i  in  I  let 

q 


j  be  the  smallest  column  index  such  that  (i,j)  «  T;  and  let  be  the  set  of 

★ 

all  such  cells.  Then  x.  “  r.,/y.  for  (i,j)  r.  C  .  As  before  let  d.  *  y,/y  . 

i  ijj  q  i%  j  q 


Removing  (p.q)  from  T  forces  X  to  0  and  we  can  calculate  the  new 

pq 


solution  x. 


and  y.  for  all  i  and  j.  From  this  we  can  find  r  as 

j  pq 


x  y  .  Once  (p,q)  is  dropped  from  T  the  primal  problem  (V)  becomes 
p  q  -  w 


Minimize 


£  x.  +  £ 

iel  1  iel 
P  9 


£  e.  xn  +  £ 

iclp  ip  P  (i.j)«Cq 


/JjLlV  ± 
\Jq  J  yq 


subject  to  the  constraint 
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(iil)  some  cell  (h,k)  in  1^  x  -  {(p,q)}  becomes  tight; 

* 

then  add  (h,k)  to  T  remove  (p,q)  from  T,  set  r  =  r  , 

r  4  r  4 

and  go  to  2. 

* 

(c)  All  cells  (p,q)eT  satisfy  X  >0  and  r  =  r  or  X  *  0 

3  pq  -  pq  pq  pq 

*  o 

and  r  =  r  >  r  .  Go  to  4. 

pq  pq  -  pq 

* 

(4)  The  optimum  tree-restricted  solution  to  R  is  given  by  the 
current  T,  x^,  y^  and  X  t  since  these  satisfy  the  duality  theorem.  The 
optimum  forest  F  for  R  is  obtained  by  dropping  those  cells  (p,q)  from 
T  with  X  =0;  the  same  x. .  y.  and  X..  are  optimal  for  R  with  forest 

pq  i  tj 


THEOREM  7.  For  a  non-degenerate  problem  the  algorithm  converges  in  an 
infinite  number  of  steps  to  the  optimum  answer  to  problem  R. 

A  proof  of  this  theorem  appears  in  Smeers  [9] . 

An  intuitive  proof  is  the  following.  Suppose  the  algorithm  always  keeps 
the  same  tree;  i.e.,  steps  2  (a)  (ii)  and  2  (b)  (iii)  are  never  entered.  Then 
it  is  easy  to  show  that  the  tree  restricted  problem  (H),  after  making  the  sub¬ 
stitutions  **  rtj/y  for  (i,j)cT,  is  a  concave  problem,  and  the  algorithm 
can  easily  be  shown  to  converge  (infinitely),  see  Zangwill  [12].  Clearly  there 
are  only  a  finite  number  of  possible  tree  structures,  and  each  one  has  (in  the 
non-degenerate  case)  a  different  optimum  value.  It  follows  that  steps  2(a)  (ii) 
and  e(b)  (iii)  of  the  algorithm  will  be  entered  only  a  finite  number  of  times  so 
that  the  processes  eventually  settles  down  on  a  single  tree  for  which  the  pre¬ 
vious  argument  holds. 


Since  the  above  algorithm  involves  finding  square  roots  at  each  step  which 
are  implemented  by  algorithms  that  converge  infinitely,  there  can  be  no  strictly 
finite  procedure  for  solving  the  above  problem.  However,  the  following  modifica¬ 
tion  of  step  (2)  of  the  algorithm  ensures  that  all  the  rest  of  the  calculations 


are  finite. 
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(2  )  Examine  Che  cells  (p.q)  in  the  current  basis  T  to  find  the  set  S  of 

all  cells  (p , q)  such  that  either  (a)  X  <  0  or  (b)  X  >0  and  r"  >  r  . 

pq  pq  pq  pq 

If  S  =  0  go  to  (4).  Otherwise  let  F  *  T  -  S  and  solve  the  resulting  forest 
restricted  problem  (H)  by  a  Lagrangian  procedure  similar  to  that  in  the  proof  of 
Theorem  5.  If  x  y  >  r  for  all  (p,q)eS  then  the  optimum  fore  st  restricted 

p  q  pq 

solution  has  been  found.  If  not,  add  cells  of  S  to  F  in  all  possible  ways 
and  resolve  until  the  condition  holds.  A  branch  and  bound  procedure  can  be  de¬ 
vised  to  simplify  this  calculation. 

It  is  easy  to  show  that  with  step  (2  ')  the  algorithm  becomes  a  finite  one 

(except  for  the  square,  roots)  since  there  are  only  a  finite  number  of  forest 

bases,  and  step  (2  )  finds  the  minimum  objective  value  for  each  one  when  it 
is  considered.  Hence  no  forest  basis  can  be  considered  more  than  once  in  the 
computation  process. 

Since  computers  have  limited  accuracy,  infinitely  convergent  processes 
have  to  be  (finitely)  terminated  after  they  converge  to  that  accuracy.  The 
authors  have  found  that  the  original  algorithm  with  step  (2)  instead  of  (2') 
converges  quickly  to  within  the  accuracy  limit  determined  by  the  computer. 

In  fact,  the  number  of  pivet  steps  needed  seems  to  be  fewer  than  that  needed 
for  a  corresponding  transportation  pi'oblem  (see  [11]).  Therefore,  step  (2') 
has  noc  (as  yet)  been  programmed. 

In  our  codo  we  have  implemented  degeneracy  prevention  techniques  similar 
to  thosi'  used  in  linear  programming, 

Smeers  [  q]  has  proposed  an  alternate  way  of  finite  termination  for 


the  algorithm. 
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We  now  discuss  two  way6  of  implementing  the  starting  routine  of  the 
algorithm.  The  first  procedure  is  similar  to  the  improvement  routine  of  the 
algorithm. 

Starting  Routine  1.  Find  an  initial  tree  basis  T  by  any  means.  A 

good  heuristic  is  to  try  and  get  as  many  of  the  large  entries  in  R  as 

possible  into  this  initial  basis,  Now  solve  for  x..  and  y  using  T. 

* 

For  each  non  basic  cell  (i,j)  if  x^^  y^  <  r  replace  r^  by  r_  =  xi^j* 

After  this  has  been  done  go  back  and  for  each  (i,j)  such  that  r  <  r^ 

i k 

increase  r.  ,  until  either  a  new  ceil  becomes  tight  and  enters  the  basis 
1 1 

k 

in  place  of  (i,j)  or  else  r. .  =>  r, ..  Note  that  this  makes  g  constantly 

ij  ij 

increase  and  hence  these  steps  are  just  the  reverse  of  the  improvement 

routine  of  the  algorithm.  After  a  finite  number  of  such  steps  a  primal 

feasible  tree-restricted  solution  to  the  original  problem  will  be  attained. 

Starting  Routine  2.  Select  an  arbitrary  set  of  positive  y/s  such 

that  I",  v  *  1.  A  good  choice  would  be  to  select 

y  -  (1-,,1/CIt  ) 

1  i  j 

Now  choose  x.  =  Max  (r. ./y.)  and  put  all  tight  cells  into  the  basis.  If 
i  .  ij  7j 

there  is  a  column,  say  q.  with  no  tight  cells,  select  an  arbitrary  row, 

say  p,  and  raise  r  to  the  value  x  y  and  add  this  cell  to  the  basis  also. 

pq  p  q 

We  now  have  a  primal  feasible  forest  basis  which  can  further  be  extended  to 
a  primal  feasible  tree  basis  by  using  the  techniques  of  the  proof  of 
Theorem  3. 

As  a  final  remark,  we  would  like  to  discuss  how  a  primal  feasible 

tree-solution  can  be  used  to  determine  the  next  tree  solution.  For  the 

first  solution  we  have  y,  *  d.,y,  for  k  fixed  and  all  t.  If  we  now 

'j  jk  k 

change  r  for  (p,q)  being  a  tight  cell,  the  values  of  d  , 
pq  jk 


will  be 
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affected  only  for  those  y  such  that  (p,q)  is  a  part  of  the  unique  path 

it  irk 

from  column  j  to  column  k.  If  r  is  changed  to  r  then  either 

pq  pq 


*  ** 

d*  =  d  -E£  or  d*  *  d  ■ 
jk  jk  r**  jk  djk  r* 

pq  pq 


depending  on  where  (p,q)  is  in  the  path  from  column  j  to  column  k. 


6 .  EXAMPLES 

Tha  first  example  is  designed  to  demonstrate  most  of  the  steps  of  the 
algorithm.  We  start  with  the  data  and  an  initial  tree  as 


Cells  (3,3)  and  (1,3)  are  not  primal  feasible  since 


Following  Starting  Routine  1  we  replace  the  problem  by 


whose  value  is  g  =  51.3.  We  now  bring  cell  (3,3)  into  the  basis  and  can 

remove  any  cell  in  Tj  -  {(2,3),  (1,2),  (3,1)}.  We  choose  to  remove  cell  (3,1) 

* 

and  now  increase  r^  ^  from  2.8,  trying  to  raise  it  to  4,  without  causing 
primal  infeasibilities.  We  succeed  and  obtain  the  problem: 
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vhose  value  is  g  *  55.93.  We  now  bring  cell  (1,3)  into  the  basis  and  can 

remove  either  cell  in  *  {(1,2),  (2,3)}  we  choose  to  remove  (2,3) 

(although  later  it  will  become  evident  that  the  other  choice  is  better)  in 

* 

order  to  demonstrate  more  steps  of  the  algorithm.  We  raise  r^  to  4 
without  encountering  primal  infeasibilities  and  obtain  the  problem 


© 

© 

© 

3 

© 

7 

4 

3 

whose  value  is  g  *  56.  It6  primal  feasible  solutions  are 

/  5  5  4  \ 

y  “  ll4  ’  14  ’  14/ 

x  -  (l4,  28,  14) 

'  /  1  JL  I) 

X  \  4*  2’  J 

and  the  dual  solutions  are  given  by 


5 

t 

l 

11  ’  yl  “  14’ 

P22 

m  v 

x2 

*  2  * 

P33 

5 

1 

2 

12  *  y2  "  p22 

"  14 

■  —  • 
2 

'  14* 

Pi  3 

Since  (1,2)  is  the  only  cell  with  negative  dual  variable  we  have  p  ■  l,  q  *  2, 
Ix  -  {2},  J2  -  {1,3},  R2  =  {(1,1), (3, 3)},  -  {(2,2)}  and  en  -  e13  -  1, 

d22  =  1.  Hence  we  can  calculate  from  Theorem  5  as 


1  1 


12 


10 

1 

1 


-  /4T-  6.71 


We  also  have  Ii  x  ^  2  *  {(2,1),  (2,3)}  as  the  two  cells  that  may  become  tight 
as  we  increase  r^.  For  them  we  have: 
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16.67  from  cell  (2,1) 


and 

*  10*4 

r12  -  ”~7  *  5.71  from  cell  (2,3). 

The  smallest  constraint  comes  from  cell  (2,3)  so  we  bring  it  into  the  basis 
obtaining  the  problem: 


© 

5 

© 

3 

© 

© 

4 

3 

© 

whose  value  is 

8  -  55. 

18.  Its 

solution  is 

/'  35 

40 

28 

\ 

y  ■ 

Is 

103  » 

103 

) 

X  » 

103 

103 

\ 

1 

\  7  » 

4  » 

7 

J 

/ 

v  m. 

7 

4 

\ 

A  * 

V  15  * 

15  * 

15 

; 

The  dual  solution  is  given  by: 

35 


J11 


yl  =  103  ’ 


y22 


40 

103 


p33  "  x3 


_4_ 

15 


P23  *  x2 


40 


P22  15  103  *  ,08,  pi3 


X1  '  P11  "  15 


35 

103 


-.o; 


Cell  (1,3)  is  the  only  one  with  negative  dual  variable  90  p  ■  1,  q  ■  3,  Ij 
J3  *  £l}»  Rj  *  {(1,1)},  -  {(2,2),  (3,3)},  and  ^  -  1,  ^  , 


d33  "  X* 


Hence  we  calculate 


10-fe+  f 


5  x 


4.76 
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We  also  have  Ij  x  ■  {(2,1),  (3,1)}  so  that  the  other  constraints  on 

* 

r^  are  given  by 


* 

5*7 

13  - 

3 

* 

4*5 

13  - 

4 

11.67 


Therefore  we  increase  r^  to  4, 


from  cell  (2,1) 
from  cell  (3,1) 

76  obtaining  the  problem 


® 

5 

(Tt?) 

3 

© 

© 

4 

3 

© 

whose  value  is  g  =»  54.83.  The  primal  and  dual  solutions  are: 
y  =  (.29,  .41,  .30) 
x  -  (16.56,  24.35,  13.92) 
x'  *  (.30,  ,44,  .26) 

P11  “  yl  *  *30’  Pi3  3  xi  "  Pn  3  °>  P22  *  y2  “  *41 

p23  *  X2  '  P22  =  ,44  '  *41  =  *03,  P31  3  x3  3  *26> 


Since  all  cells  satisfy  3(c)  of  the  algorithm,  we  have  the  optimum  solution  to 

•jjf 

both  R  and  R  .  We  next  solve  Evans'  problem  [4],  which  is  given  with  an 
initial  tree  basis: 


15 

© 

© 

© 

13 

0 

15 

17 

© 

© 

-J2-J 

J2 _ 

cells  (2,2)  and  (4,2)  are  not  primal  feasible  since 


13  > 


13  *  22  •  23 


34  •  44 


«  4.4 


We  now  bring  cell  (2,2)  into  the  basis  and  remove  cell  (2,1).  The  problem 
remains  primal  feasible  and  we  can  increase  the  value  of  r 22  to  13  while 
remaining  primal  feasible. 


(4,3).  We  then  increase  r^2  to  12  remaining  primal  feasible. 
Initial  primal  feasible  solution  to  the  original  problem: 


x  -  (132.17,  74.70,  105.13,  68.96) 
x'  -  (.35,  .20,  .28,  .18) 
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Dual  variables 

P41  "  yl  =  -49 


P42  =  X4  *  P41  =  -18  ‘  ‘49  =  ‘-31 


p33  x3  *28 


P22  =  *2  =  '2° 


13 


*  y3'  -  p33  -  .33  -  .28  =  .05 


pl7  ~  X1  "  p13  *  =  *3° 

Only  cell  (4,2)  has  a  negative  dual  variable. 

Take  cell  (4,2):  I  x  J  =  [(1,1),  (2,1) 

9  P 

For  (1,1)  15  <  23  •  34  so  r*  <52.1 


42 


For  (2,1)  13  <  13  •  34  so  r*2  <  36.8 


42 


For  (3,1)  15  <  35  •  23  •  34 


so  r 


44 


42 


42 


Therefore  r^  <  36.8  but  r°2  =  25.18 


Therefore  we  can  raise  r . „  to  25.18 

42 


15 

© 

© 

13 

© 

0 

15 

17 

r© 

© 

(STTe) 

22 

(.32,  .23,  .45) 

(98.00,  55.43,  78.01,  107.33) 
(.29,  .16,  .23,  .32) 


(3,1)} 


41.4 
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Dual  variables: 

P41 

“  yl 

=  .32 

P42 

/ 

-  X, 

4 

-  P41  = 

/ 

x3 

-  .23 

p22 

/ 

=  X2 

=  .16 

p13 

=  y3 

-  P33  -  - 

.23  =  .22 

p12 

/ 

=  X1 

p13  =  ’ 2 9 

.22  =  .07 

Since  we  satisfy  3(c)  of  the  algorithm,  we  havr  the  optimal  :oiuti-n 

-JL 

to  R  and  R  . 

7.  CONCLUSIONS 

The  algorithm  presented  in  this  paper  has  several  advantage .>  . 

a. ,  *£  we  most  stop  before  the  optimum,  solution  to  the  problem  i  •, 

reaci.  d,  the  non- line "rity  makes  g  Close  to  the  optimum.  Thi. 
solution  will  be  much  closer  tnan  would  be  the  ca  ■■■<■  in  a  corres¬ 
ponding  similar  linear  problem. 

b. )  The  method  always  keeps  a  primal  ‘easi.h’.c  solution  so  tiiat  one 

can  stop  the  procedure  at  any  tins:  and  havi  a  usable  solution. 

c. )  In  order  to  find  the  negative  restricted  dual  variable,  only  the 

basis  tree  must  be  searched.  A  most  negative  indicator  rule 
would,  therefore,  be  available  at  low  computational  cos  .  (This 
is  not  the  case  tor  transportati on  problem  (see  [111), 

d .  )  The  search  icr  new  limiting  tells  requires  searching  on  ch 

areas  1  x  J  or  1  x  J  ,  anH  not  the  whole  matrix. 

p  q  q  r 

c.)  Previous  solutions  to  the  jrimal  problem  can  be  u'-ed  to  enerate 


feasible  solutions  co  probl 


wi th  simi lar  data. 
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£.)  Accuracy  is  not  a  problen  since  a  solution  to  a  forest  basis 
can  be  found  independently  of  any  previous  op  ;r;tions. 

Retention  of  a  tree  solution  throughout  the  computation  and  use  of 
previous  solutions  makes  the  steps  in  the  algorithm  very  similar  to  those  if 
transportation  problems.  Srinivasan  and  Thompson  report  excellent  computa¬ 
tional  results  of  175  x  175  transportation  problems  in  sc”en  seconds  [11]. 
The  number  of  pivots  required  for  the  modular  design  problems  tested  so  fa. 
are  somewhat  fewer  than  for  a  transportation  problem  of  the  same  size.  This 
means  that  modular  design  problems  can  be  solved  in  only  slightly  more 
computer  time  than  that  required  for  comparable  transportation  pr  '"ilenia . 

The  authors  are  currently  preparing  a  report  on  computational  experience 
with  their  modular  design  .-ode. 
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